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Abstract 

We investigate the local regularity of pointed spacetimes, that is, time- 
oriented Lorentzian manifolds in which a point and a future-oriented, 
unit timelike vector (an observer) are selected. Our main result covers 
the class of Einstein vacuum spacetimes. Under curvature and injectivity 
bounds only, we establish the existence of a local coordinate chart defined 
in a ball with definite size in which the metric coefficients have optimal 
regularity. The proof is based on quantitative estimates, on one hand, for 
a constant mean curvature (CMC) foliation by spacelike hypersurfaces 
defined locally near the observer and, on the other hand, for the metric 
in local coordinates that are spatially harmonic in each CMC slice. The 
results and techniques in this paper should be useful in the context of 
general relativity for investigating the long-time behavior of solutions to 
the Einstein equations. 

1 Introduction 

1.1 Quantitative estimates for CMC foliations 

We denote by (M, g) a spacetime of general relativity, that is, a time-oriented, 
{n + l)-dimensional Lorentzian manifold whose metric g, by definition, has 
signature (— , . . . , -f ). Our main result in the present paper will concern 
vacuum spacetimes, that is, Ricci-flat manifolds, although this assumption will 
be made only later in the discussion. Building on our earlier work [12] . we 
continue the investigation of the local geometry of Einstein spacetimes, using 
here techniques for partial differential equations. Our main objective will be, 
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under natural geometric bounds on the curvature and the injectivity radius 
only, to establish the existence of local coordinate charts in which the metric 
coefficients have optimal regularity, that is, belong to the Sobolev space 
for all real a G (1, cx)). The construction proposed in the present paper is local in 
the neighborhood of a given "observer" and, in turn, our result provides a sharp 
control of the local geometry of the spacetime at every point. This optimal 
regularity theory should be useful for tackling the global regularity issue for 
Einstein spacetimes and investigating the long-time behavior of solutions to the 
Einstein equations. 

As in [12], we consider a pointed Lorentzian manifold (M, g, p, Tp), that 
is, a time-oriented Lorentzian manifold supplemented with a point p € M and 
a future-oriented timelike vector Tp at that point. The pair (p,Tp) is called 
a (local) observer and is required for stating our curvature and injectivity ra- 
dius bounds for some given curvature constant A and injectivity radius constant 
A > 0; see (|1.2|) in Section [27T] below. We will establish the existence of a neigh- 
borhood of the observer (p, Tp) whose size depends on A, A only and in which 
local coordinates exist in such a way that the regularity of the metric coefficients 
can be controlled by the same constants. Since the Riemann curvature involves 
up to two derivatives of the metric it is natural to search for an estimate of the 
metric in the W^ "^ norm, and this is precisely what we achieve in the present 
paper. 

We will proceed as follows. Our first task is constructing a constant mean 
curvature (CMC) foliation by spacelike hypersurfaces, which is locally defined 
near the observer and satisfies quantitative bounds involving the constants A, A, 
only; see Theorem 12.21 below. Our method can be viewed as a refinement of 
earlier works by Bartnik and Simon (covering hypersurfaces in Minkowski 
space) and Gerhardt [131 [TS] (global foliations of Lorentzian manifolds) . If one 
would assume that the metric g admits bounded covariant derivatives of suffi- 
ciently high order of regularity, then the techniques in [9l [14] would provide the 
existence of the CMC foliation and certain estimates. Hence, the construction 
of a CMC foliation on a sulEciently smooth manifold is standard at small scales. 

In contrast, in the framework of the present paper only limited differentiabil- 
ity of the metric should be used and uniform bounds involving the curvature and 
injectivity radius bounds, only, be sought. We have to solve a boundary value 
problem for the prescribed mean curvature equation in a Lorentzian background 
and to establish that a CMC foliation exists in a neighborhood (of the observer) 
with definite size and to control the geometry of these slices in terms of A, A, 
only. A technical difficulty in this analysis is ensuring that each hypersurface 
of the foliation is uniformly spacelike and can not approach a null hypersurface. 
Deriving a gradient estimate for prescribed curvature equations requires the use 
of barrier functions determined from (parts of) suitably constructed geodesic 
spheres. 
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1.2 Earlier works 



An extensive study of (sufficiently regular) spacetimes admitting global folia- 
tions by spatially compact hypersurfaces with constant mean curvature is avail- 
able in the literature. Andersson and Moncrief [3 [5] and Andersson [5J [5] have 
established global existence theorems for sufficiently small perturbations of a 
large family of spacetimes. For instance, their method allowed them to estab- 
lish a global existence theorem for sufficiently small perturbation of Friedmann- 
Robertson- Walker type spacetimes. Their construction is based on construct- 
ing a global CMC foliation and uses harmonic coordinates on each slice. In 
these works, the authors derive (and strongly rely on) a priori estimates which 
are based on the so-called Bel-Robinson tensor and involve up to third-order 
derivatives of the metric. In contrast, we focus in the present paper on the local 
existence of such foliations but require only the sup norm of the curvature to be 
bounded. Our new approach leads to a construction of "good" local coordinates 
(see below) and allows us to explore the local optimal regularity of Lorentzian 
metrics. Another direction of research on CMC foliations is currently developed 
by Reiris [26l |27] , who analyzes the CMC Einstein flow in connection with the 
Bel-Robinson energy and also imposes higher regularity of the metric. 

We also refer the reader to an ambitious program (the L"^ curvature con- 
jecture) initiated and developed by Klainerman and Rodnianski in a series of 
papers; see [HI [501 IHl HI] ■ In these works, the authors are interested in con- 
troling the geometry of null cones which may become singular due to caustic 
formation. The regularity of null cones is needed in order to suitably extend 
the methods of harmonic analysis to the Einstein equations. In particular, the 
recent result [22] provides a breakdown criterion for solutions to the Einstein 
equations. In comparison with the present work, the objectives in |22| are differ- 
ent: these authors rely on hyperbolic techniques and investigate the geometry of 
light cones, while our approach in the present paper is purely elliptic in nature 
and addresses the geometry of the spacetime itself. 

1.3 CMC— harmonic coordinates of an observer 

Our second task is constructing local coordinates. In Riemannian geometry it is 
well-known that geodesic-based coordinates and distance-based coordinates fail 
to achieve the optimal regularity of the metric. The use of harmonic coordinates 
on Riemannian manifolds was first advocated by De Turck and Kazdan 13J and, 
later, a quantitative bound on the harmonic radius at a point was derived by 
Jost and Karcher |18j in terms of curvature and volume bounds, only. More 
recently, the issue of the optimal regularity of Lorentzian metrics was tackled 
by Anderson in the pioneering work [51 H] . He proposed to use a combination 
of normal coordinates (based on geodesies) and spatially harmonic coordinates, 
and derived several uniform estimates for the metric coefficients. This construc- 
tion based on geodesies does not lead to the desired optimal regularity, however. 
We also refer to earlier work by Anderson T, 'S] for further regularity results 
within the class of static and, more generally, stationary spacetimes. 
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Our main result covers Einstein vacuum spacetimes, that is, manifolds sat- 
isfying the Ricci-flat condition 

Ricg = 0, (1.1) 

and the construction we propose is as follows. Relying on our quantitative 
estimates for CMC foliations near a given observer ('Theorem l2.2p and then ap- 
plying Jost and Karcher's theorem for Riemannian manifolds 18J, we construct 
(spatially) harmonic coordinates on each spacelike CMC slice. We refer to such 
coordinates as CMC-harmonic coordinates, and we prove first that, on every 
slice, the spatial metric coefficients gy belong to the Sobolev space and 
satisfy the quantitative estimate 

l|g»jllp,Tp,w2,» < C(a,A, A) 

for all a < GO and some constant C(a, A, A) > (depending also on the dimen- 
sion n). In addition, we also control the lapse function and the shift vector 
associated with these local coordinates. The shift vector, denoted below by ^, 
arises since coordinates are not simply transported from one slice to another 
but are chosen to be harmonic on each slice. The lapse function, denoted below 
by A, is a measure of the distance between two nearby slices. 

In turn, we arrive at the following main result of the present paper. 

Theorem 1.1 (CMC-harmonic coordinates of an observer). There exist con- 
stants < c{n) < c{n) < 1 and C{n),Cq{n) > depending upon the dimen- 
sion n (and some exponent q £ [l,oo)) such that the following properties hold. 
Let (M, g, p, Tp) be an {n -\- I) -dimensional, pointed, Einstein vacuum space- 
time satisfying the following curvature and injectivity radius bounds at the scale 
r > 0: 

Rl,^^{M, g, p, Tp) < r-2, Inj(M, g, p, Tp) > r. (1.2) 

Then, there exists a local coordinate system x — {t,x^, . . . ,x") having p = 
(ri, 0, . . . , 0) for some ri £ [c{n)r, c{n)r] and defined for all 

\t - ri\ < c{n)\, {{x^)^ + . . . + {x")^)'^^ < c{n)\, 
so that the following two properties hold: 

i) Each slice = {{x^)^ -I- ... -1- (x")-^ < c(n)'*r^} on which t remains con- 
stant is a spacelike hypersurface with constant mean curvature c{n)~^r~'^t 
and the coordinates x :— (x^, . . . , x") are harmonic for the metric induced 
on Et. 

ii) The Lorentzian metric in the spacetime coordinates x = (t,x^ , . . . , x") has 
the form 

g = -X{xf {dtf + gij{x){dx' + C{x)dt){dx' +S,^{x)dt) (1.3) 
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and is close to the Minkowski metric in these local coordinates, in the 
sense that 

and for each q £ [1, oo) 

The theorem above estabhshes the existence of locally defined CMC-harmonic 
coordinates near any observer. The coordinates cover a neighborhood of the base 
point, whose size is of order r in the timehke and in the spacelike directions. 
In the statement above, dx^ and dxx^ denote any spacetime first- and second- 
order derivatives of the metric coefficients in the local coordinates, respectively, 
while dwst denotes the volume form induced on Et by the spacetime metric and 
can be computed in terms of the spatial coordinates x. 

Finally, let us put our results in a larger perspective. The proposed frame- 
work relies on constructing purely local CMC-harmonic coordinates and, there- 
fore, applies to spacetimes which need not admit a global CMC foliation. Hw- 
ever, based on our local regularity theory we can also control the global geometry 
of the spacetime, as follows. 

Given a pointed Lorentzian manifold (M, g, p, Tp) satisfying a global version 
of the curvature and injectivity radius estimates ()1.2p . we can find a global atlas 
of local charts covering the whole of M and in which the metric coefficients have 
the optimal regularity. Such a conclusion is achieved by introducing a notion of 
global CMC-harmonic radius viewed by the observer (p, Tp): it is the "largest" 
radius r > such that tocai CMC-harmonic coordinates exist in a ball of radius 
r about each point (and satisfy the uniform estimates stated in Theorem 11.11 
for some fixed constants c(n), c(n), C{n), Cq{n)). By establishing a lower bound 
on the radius of balls in which local CMC-harmonic coordinates exist at every 
point, we obtain the desired global optimal regularity. Again, this is a purely 
geometric result that involves the curvature and injectivity radius bounds, only. 
This development is a work in progress. 

Throughout this paper, we use the notation C, C", Ci, . . . for constants that 
only depend on the dimension n and may change at each occurrence. 

2 CMC foliation of an observer 
2.1 Main statement in this section 

In this section, we derive quantitative bounds on local CMC foliations for a gen- 
eral class of Lorentzian manifolds which need not satisfy the Einstein equations. 
For background on Riemannian or Lorentzian geometry we refer to [lOl [17l [23] . 
Let (M, g) be a time-oriented, (n + l)-dimensional Lorentzian manifold, and 
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let V be the Levi-Civita connection associated with g. Given a point p G M, 
we want to construct a constant mean curvature foUation that is defined near 
p and whose geometry is uniformly controled in terms of the curvature and 
injectivity radius, only. The inner product of two vectors X, Y is also written 
{X,Y), = {X,Y)=giX,Y). 

In fact, rather than a single point on the manifold we must prescribe an 
observer, that is, a pair (p,Tp) where Tp is a unit, future-oriented, timelike 
vector at p (also called a reference vector). We use the notation (p, Tp) £ 
Tj^M for the bundle of such pairs, and we refer to (M, g, p, Tp) as a pointed 
Lorentzian manifold. The vector Tp naturally induces a (positive-definite) inner 
product on the tangent space at p, which we denote by gTp = ( • , • )tp- We 
sometimes write |^|gr for the Riemannian norm of a vector X. To simplify the 
notation, we often write T instead of Tp . 

On a Lorentzian manifold the notion of injectivity radius is defined as follows. 
Consider the exponential map expp at the point p, as a map defined on the 
Riemannian ball B^.^ (p, r) (a subset of the tangent space at p) and taking values 
in M; this map is well-defined for all sufficiently small radius r, at least. 

Definition 2.1. Given a Lorentzian manifold (M, g), the injectivity radius 
Inj(M, g, p, Tp) of an observer (p, Tp) G T^M. is the supremum over all radii 
r > such that the exponential map expp is well-defined and is a global diffeo- 
morphism from the subset B^.^. (p, r) of the tangent space at p to a neighborhood 
of 'p in the manifold denoted by Sg.j,(p,r) :— expp(i?g.j,(p, r)) C M. 

To simplify the notation, we will also use the notation i?T(p, r) and St(p, r) 
for the above Riemannian balls. To state our assumption on the curvature we 
need a Riemaniann metric defined in a neighborhood of the point p. This 
reference metric is also denoted by gx and is defined as follows. 

By parallel transporting the vector Tp, with respect to the Lorentzian con- 
nection V and along radial geodesies leaving from p, we construct a vector field 
T which, however, may be multi-valued since two distinct geodesies leaving 
from p, in general, may eventually intersect. We use the notation T^ for the 
vector field defined along a radial geodesic 7 lying in the set St(p, f")- Then, 
to this vector field we canonically associate a positive-definite, inner product 
gT^ = ( • I • )t^ defined in the tangent space of each point along the geodesic. 
We write \A\t or \A\t for the corresponding Riemannian norm of a tensor A. 

We then consider the Riemann curvature Rm of the connection V and, 
given an observer (p, Tp), we compute its norm in the ball of radius r 

RLax(M, g, p, Tp) := sup |Rm|T, , (2.1) 

7 

where the supremum is taken over every radial geodesic from p of Riemannian 
length r, at most. Note that |Rm|T^ is evaluated with the Riemannian reference 
metric rather than from the Lorentzian metric. Our main assumption (jl.2p is 
now well-defined. 

Our objective in this section is constructing a foliation near p, say (J^^j^j St, 
by n-dimensional spacelike hypersurfaces St C M , and we require that each 
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slice has constant mean curvature equal to t, the range of t being specified by 
some functions t = i(p) and t = i(p). Moreover, this foliation should cover a 
"relatively large" part of the ball Bt{p, r). 

Theorem 2.2 (Uniform estimates for a local CMC foliation of an observer). 

There exist constants c, c, c, 0, G (0, 1) with c < c <c, depending only on the di- 
mension of the manifold such that the following property holds. Let (M, g, p, Tp) 
be a pointed Lorentzian manifold satisfying the curvature and injectivity radius 
assumptions (jl.2[) at some scale r > 0. Then, the Riemannian hall ST(p,cr) 
can be covered by a foliation by spacelike hypersurfaces with constant mean 
curvature t, 

( U St) DST(p,cr), t:^n^—^, < ;= „i±i, ^2.2) 

t<t<t 

in which the time variable describes a range [t, t\ determined by some real s € 
[c, c] and, moreover, the unit normal vector N and the second fundamental form 
h of the foliation satisfy 

1 < -5(N,T) < 1 + 6l-\ 9 < -r-^g{Vt,Vt) < e-\ r\h\<e-\ 

(Recall that the vector field T is defined by parallel translating the given vector 
Tp along radial geodesies from p.) 

Hence, a foliation exists in a neighborhood of the base point, in which the 
time variable is of order 1/r and describes an interval with definite size. In our 
construction given below, it will be important that s be chosen to be sufficiently 
small. 

Note that the above theorem is purely geometric and does not depend explic- 
itly on the coordinates that we are going now to introduce in order to establish 
the existence of the above foliation and control its geometry. 

The rest of this section is devoted to giving a proof of Theorem l2.2l We first 
construct the CMC hypersurfaces as graphs over geodesic spheres associated 
with the Lorentzian metric. Geodesic spheres associated with the reference 
Riemannian metric will be introduced to serve as barrier functions. Indeed, 
each CMC hypersurface will be pinched between a Lorentzian and a Riemannian 
geodesic ball. The level set function describing the CMC hypersurface satisfies 
a nonlinear elliptic equation, whose coefficients have rather limited regularity, 
and this will force us to use the Nash-Moser iteration technique. 

2.2 Formulation in normal coordinates 
Foliation by geodesic spheres 

We begin by introducing spacctimc normal coordinates and by expressing the 
prescribed mean curvature equation in these coordinates. As we established 
earlier in [121, under the curvature and injectivity radius assumption (|1.2p for the 
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observer (p, Tp), there exist positive constants c < c <c < 1 and C depending 
only on the dimension of the manifold, such that the following properties hold. 

First of all, the foliation by subsets Tir of geodesic spheres is defined as 
follows. Let 7 be a future-oriented, timelike geodesic containing p and let us 
parameterize it so that p = j{cr). Set q := 7(0) and consider the (new) observer 
(q, Tq) with Tq := 7'(0). The constant c is chosen sufficiently small so that the 
injectivity radius of the map expq (computed for the observer (q, Tq)) is cr, 
at least. From the point q we consider normal coordinates y = (y") = {T,y^) 
determined by the family of future-oriented timelike radial geodesies from q, 
so that the Lorentzian metric takes the form g — —dr^ + gijdy^dy^. These 
coordinates cover a part of the future of the point q and at least the region 

C+(q,cr) := expq(C+(q,cr)), 

C+(q,cr) := |yei?TjO,cr), gT,(y,V) < 0, ^^^^^^yy], > ^ -c^ ■ 

The base point p is identified with (r, . . . , y") = (cr, 0, . . . , 0) in these co- 
ordinates. By relying on the curvature bound, analyzing the behavior of Jacobi 
fields, and using standard comparison arguments from Ricmannian geometry, 
one can establish in well-chosen coordinates 1121: 



C ^Sij < < C Stj, 



9gj, 



dr 



■V 



(2.3) 

<C in C+(q,cr) n {t < r < r}. 



where r := cr and r :— cr. 

The reference Riemannian metric associated with the vector field d/dr (ob- 
tained by parallel transporting the vector Tq) reads g := dr^ -I- gij dy^dy^ . We 
use the notation d(-,-) and S{-,-) for the distance function and the geodesic 
spheres associated with the metric g, respectively. By construction, the func- 
tion T coincides with the distance function d(q, •). It will be useful also to have 
the following estimate of the Riemann curvature of the Lorentzian metric 

|Rm|g <Cr-2 in C+(q, cr) n {r< r <r}, (2.4) 

which is based on the reference metric g. 

In turn, the above construction provides us with a foliation (by Lorentzian 
geodesic spheres) of some neighborhood of p (with definite size) by n-dimensional 
spacelike hypersurfaces Ht, hence p e |J^gj^_j Ht- 

Now, consider the time fmiction t. The standard Hessian comparison theo- 
rem for distance functions in Ricmannian geometry is also useful in Lorentzian 
geometry and, more precisely, shows that the (restriction of the) Hessian of r is 
equivalent to the induced metric: 

k{T,r)g,, < (-VV)|£,y <'fc(T,r)gy, (2.5) 
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where E := (Vr)"*" denotes the orthogonal complement and 



r-WC - r-WC 

k{T,r) := , k(T,r) 



tan(Tr-i\/C) ' t&nh {t r-^VC)' 

Observe that both fc(T, r) and k{T,r) behave like 1/t when r — s- 0. Note also 
that k will remain non-singular within the range of interest, since t will be 
chosen to be a small multiple of r. 
Consequently, noting that 

-nr.i^=:^„ (2.6) 

taking the trace in the inequalities (j2.5p . and then using the uniform estimate 
(|2.3p . we see that the mean curvature of each slice r = const, is close to n/r. 
Our objective now is to replace these slices by constant mean curvature slices 
by making a small perturbation determined by solving an elliptic equation in 
these normal coordinates. 

We will also use the Riemannian geodesic spheres associated with g. Namely, 
consider the Riemannian distance function d(p', •) computed from some arbi- 
trary point p' := 7(t) with t G [z,t]. From the expression of the reference 
metric we find _ 

V^d(p',.) = Vd(p',.)-2|^A (2.7) 

where V is the covariant derivative associated with g. Again by the Hessian 
comparison theorem and since |9d/9T| < |Vd|g = 1, we find after setting 

E := (Vd)^ 

Md,r)-^) g\g<{V'd)\s<(k(d,r) + ^^ g|^. (2.8) 

Choosing now the time variable to be a (small) multiple of r and taking 
the trace of the above inequalities, we deduce that for any a e [cr, cr] the mean 
curvature if^(p',a) (computed with respect to the ambient Lorentzian metric) of 
the (future-oriented, spacelike, and possibly empty) intersection A(jp',a) of the 
g-geodesic sphere S{p',a) and the future set C+(q, cr) satisfies the inequalities 

nk{a,r) < Hji^(^pi < nk{a,r), a<E[cr,cr\. (2-9) 

Hence, the mean curvature of the Riemannian slices enjoys the same inequalities 
as the ones of the Lorentzian slices TLt- Later in this section, we will use the 
graph of the Riemannian geodesic spheres as barrier functions. 

This completes the discussion of a domain of coordinates y covering a neigh- 
borhood of p, in which we can assume that all of the above estimates are valid. 
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Mecin curvature operator 

We now search for a new foliation IJ^ Et in which the hypersurfaces have constant 
mean curvature and can be viewed as graphs, say := {G*{y) := , 
over a geodesic leaf Ht for a given value r of time-function. Here, t is a real 
parameter varying in some interval of definite size and the functions y u*{y) 
need to be determined. In the following, we often write E = Ej, u = u*, 
and G = G*. Setting uj := du/dy^, the induced metric and its inverse are 
determined by projection on the slice E and read 



9ij =Sij -UiUj, g 



1- |Vu|2 ' 

and the hypersurface S is Riemannian if and only if 

\Vuf = g'^{u,-)uiUj < 1. 

We are interested here in spacelike hypersurfaces, and we denote by V the 
covariant derivative associated with the induced Riemannian metric gij. We 
easily obtain 

V«|2 ,2 |V«|2 



1_|Vm|2' I l + |Vu|2- 
The future-oriented unit normal to each hypersurface takes the form 



N = -V1 + |Vm|2(1,Vw). 

The second fundamental form of the slice E is determined by push forward (with 
the map G) of the coordinate vector fields Yj := d/dy^: 

hij := {Vg,y,G.Yj,-N) 

1 ^ 1 1 dgu 1 kt dgij 
' V,V,« + - — 1 - -g -j-UkUj - -g -jr^uuu, 



v/l-|Vw|2 V ^ ^ 2 dT 2* dT ^ ' 2^ dT 
1==(V.V,. + A.,), 



\/l+|V«|2 

(2.10) 

where 



and r*j arc the Christoffcl symbols of g. The tensor field 'Vi'Vju is the space- 
time Hessian of the function u (restricted to the hypersurface t = u), while 
ViVjU is the fully spatial Hessian defined from the intrinsic metric gij. 

The mean curvature of a slice is the trace of hij, that is, in intrinsic form 
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where A is the Laplace operator in the hypersurface, or equivalently in local 
coordinates 

Mu= r-^^ 4l ( Vg(^v)KVu)g'^(w,-)|^) 

+ (^iy{Vur'g^Hu, .) + iy{Vu) s'H^, ■)g''{u, ■)ukUi)^^{u, ■), 

where we have introduced the nonlinear function 

u{Vu) := , \^ ,^ = Vl + |V«P = uiWu). 
VI - I Vm|^ 

Note that, in fact, i/ depends also on u. 

Local formulation of the prescribed mean curvature problem 

We are now ready to introduce a formulation of the problem of interest, in terms 
of the reference Riemannian metric g. Recall that 7 is a fixed, future-oriented, 
timelike curve passing through p. Assuming for definiteness that 2c + 4c^ < c, 
from now on we fix some s € [c, 2c] and we introduce the point ps := 7((s + s^)r) 
which lies in the future of the base point p since 7 is a timelike future-oriented 
curve passing through p for the parameter value cr. We then introduce the 
subset fig C {t = sr} whose boundary is defined by the condition 

dfls := A{ps, (s^ + s^)r) D {t = sr} 

and which fills up its interior. This choice is essential for the mean curvature 
equation (discussed below) to admit the Riemannian slices as barrier functions. 
Observe that 

Bsr {l{sr), c 0« C B,r {l{sr), 2s^/^r) . (2.11) 

Here, Bsr{l{sr), a) is the geodesic ball of radius a which lies in the slice t = sr 
and is determined by the metric gjj induced on the geodesic leaf. 

Finally, given a G (0, 1) and a bounded function H of class defined on 
fls and satisfying the restriction 

nk{sr,r) < H < nk{2s^r,r), 

we seek for a spacclike hypersurface with mean curvature H and boundary dO-g. 
Analytically, this is equivalent to solving the Dirichlet problem 

Mu = H in Vis, 

■ (2.12) 
u = sr m ollg, 

in which, therefore, wc have prescribed both the boundary of the unknown 
hypersurface and its mean curvature. In the present paper, we are mainly 
interested in the case that H is a. constant function. We also assume that s is 
sufficiently small. 
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2.3 Statements of the uniform estimates 

To establish the existence of CMC hypersurfaces as graphs over a given geodesic 
leaf 'Hr^ the main difRculty is to bound |Vm| away from 1, for all functions u 
satisfying and u = r on the boundary and having their mean curvature pinched 
in some interval. Precisely, the rest of this section is devoted to the proof of the 
following result. 

Recall that (M, g, p, Tp) denotes a pointed Lorentzian manifold satisfying 
the curvature and injectivity radius assumptions (|1.2p at some scale r > and 
that 7 is a future-oriented timelike geodesic satisfying 7(cr) — p. 

Proposition 2.3 (Uniform estimates for CMC hypersurfaces). There exist con- 
stants c,9 > depending on the dimension n only such that the following prop- 
erty holds with the notation introduced in this section. For any s G [c, 2c] and 
t e [nk{sr,r),nk(2s'^r,r)] there exists a solution u to the Dirichlet problem 
(|2.12|) associated with the (constant) mean curvature function H = t and such 
that 

sup|Vu|<l-0, supr|/i| < 6I~\ 

where ft'^ = Bsr{j{sr), C {t = sr}. 

Observe that the bound on the second fundamental form holds only in a 
subset of ils, whose diameter, however, is also of the order r. Theorem 12.21 
is immediate once we establish Proposition 12.31 The proof of Proposition 12.31 
will follow from several preliminary results. The first lemma below is a direct 
consequence of the maximum principle for elliptic operators. The other lemmas 
will be established in Subsection 12.41 

Lemma 2.4 (Comparison principle). Given two functions u, w satisfying Mu > 
Mw in their domain of definition and u < w along the boundary, one has either 
u < w in the interior of their domain of definition or else u = w. In particular, 
if M.U > nk{f,r) everywhere and u < T along the boundary, then u < T. 
Similarly, if Mu < nk{T,r) everywhere and u > t along the boundary, then 
u> T. 

Lemma 2.5 (Boundary gradient estimate). For any solution u of (|2.12p with 
mean curvature function satisfying nk{sr,r) < Mu < nk{2s'^r,r) one has 

|Vm| < — on the boundary dflg- (2-13) 

Lemma 2.6 (Global gradient estimate). Under the assumptions of Lemma \2.5\ 
one has 

sup|i/(Vu)| < Ci(n), 
where the constant Ci{n) depends on the dimension, only. 
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Now, in view of Lemmas l2.5l and l2.6l and by standard arguments [TB], one can 
check that for each t G [k{sr,r),k(2s'^r,r)] the Dirichlet problem (|2.12p admits 
a smooth solution u determining a slice with constant mean curvature t. 

Note that, by Lemma 12.61 the induced metric on Et is equivalent to the 
metric gy on the domain so that we can use, for instance, Sobolev inequalities 
on St. 

Lemma 2.7 (Interior estimates for the second fundamental form). Under the 
assumptions of Lemma \2.5[ for all q G [l,oo) there exist positive constants C2{n) 
and C^{n^ q) such that for every p' G E \ 9E 

W)\< ^^(^^^ 



d(p',9E)' 

1 \ C* ( ) 



where dus is the induced volume form on T, and d(j)' , dTi) is the distance to the 
boundary dYi associated with the induced metric gij on E. 

Observe that the upper bound in the above lemma blows-up if the point 
p' approaches the boundary of the CMC slice, and that for p' G Jl^ the factor 
d{p' , 9E) is of order r, as required for Proposition [ 



Lemma 2.8 (Time-derivative of the level function). Under the assumptions of 
Lemma \2.5\ there exist constants C^ln), C^ln) > such that 

Ciiny < < C5{ny on n'^. 

2.4 Derivation of the uniform estimates 

Proof of Lemma \2.5\ We use here the maximum principle stated in Lemma l2.4l 
The part of the Riemannian geodesic sphere S'(ps, (s^ + s^)r) (defined by g) 
"below" 17s, that is the part corresponding to r < sr, is the graph y ^ {u{y), y) 
of a function u over f2s whose boundary values are sr on d^s- Since s is 
sufficiently small, one easily checks that, for instance, 

|Vu| < s^'^ on VLs (2.14) 

and, in particular, u satisfies (|2.13l) along the boundary d^s- 

Suppose now that there exists a spacelike hypersurface {u{y),y) defined 
over having the same boundary values as the function u and such that its 
mean curvature H remains bounded in the interval [nfc(sr, r), nk[[s'^ +s^/'^)r^ r)]. 
Let us set 

TO= sup u{y), TO= sup d((u(y),2/),ps) 
and use the following comparison technique. 
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Note that the range of the function u lies between (s — 2s^)r and (s + 2s^)r, 
since u is spacehke (so the norm of its gradient can not exceed 1), its boundary 
value is sr, and the diameter of the set fig is As^r at most. We are going to 
show the pinching property 

u — sr < u ~ sr < 0, 

which immediately implies the desired boundary gradient estimate (I2.13p . 

First of all, we claim that m = sr. If this were not true, then the maximum 
of u would be achieved at some point ya in the interior of fls . Since the graph 
of u is below the graph of t = m and both graphs are tangent at the point 
(yo^fn), we conclude that at the point (j/g, uivo)) the mean curvature of u is less 
or equal to that of r = m. However, in view of the Hessian estimate (|2.5p this 
is a contradiction if m > sr. 

Considering next the lower bound for u, we claim that u > u on ilg. Other- 
wise, by contradiction there would exist a point yi e Qs such that (at least) 

d((u(yi), yi), Ps) = m< {s^ + 2s^)r. 

By comparing, at the base point (u(yi), yi), the mean curvature of the graph u 
and the one of the sphere 5'(q', d((u(yi), j/i), p^)) , we find that 

Mu{yi) > nk{{s^ + 2s^)r, r), 
which contradicts our assumption Mu < nk{2s^r,r). □ 



Proof of Lemma \2.6[ Step 1. We will first show that, for some sufficiently large 
p, the sup norm of i'(Vu) is bounded by its norm. By scaling, we may 
assume r = 1 from now on. Here, as in the rest of this paper, the main difficulty 
is making sure that all constants arising in the following arguments depend on 
the injectivity radius and curvature bounds, only. It will be convenient to work 
with the intrinsic form of the mean-curvature operator M, but the expression 
in coordinates will be also used in the end of the argument in order to control 
certain Sobolev constants. 

Recall that, on the hypersurface E, 

/^u + A]=v{yu)H, (2.15) 

where A denotes the Laplace operator on the slice and H is the prescribed 
mean curvature function. Observe that since the second fundamental form of 
the geodesic sphere is bounded, we have 



141 = 

hence 



^ dr 



\Au\<C\v{Vu)\^. (2.16) 
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Observe that the coefficients of the Laplace operator on S are nothing but metric 
coefficients on which, at this stage of the analysis, we have an L°° control, only. 

We are going to use the classical Weitzenbock identity applied to the level 
function u 

A\Vuf = 2 \V^uf + 2 (Vu, VAu) + 2 Ric(Vu, Vu). (2.17) 

First, recalling that G = (u, y) we estimate the Ricci curvature term by relying 
on the Gauss formula 

Rijki —'Rap-fsG°'GjGlGf — (hikhji — huhjk) 

+ UiUiRojko + UjUk'R.iooi + UjUiRioko — (hikhji — huhjk), 

where Gi — G^{-£-i) = ^ + The spacetime curvature being uniformly 

bounded, we find 

Kaff^sG^G^GlGfg^' > ~G {1 + \Wu\^) ga- 

= -C (1 + \Wu\^) {g,k + u^Uk). 

Taking the trace of the Gauss formula, we obtain a lower bound for the Ricci 
curvature of the hypersurface: 

Rtk > hiihkjg'^ - Hh^k - C{1 + \'^uf){gik + u^Uk), 
and therefore 

Ric(Vu, Vu) > -C(l + IVmH^ 
In turn, from (j2.17p we deduce the key inequality 

A|Vw|2 - 2|v2u|2 > 2 (Vw, V (Au)) - C {l + |Vup)^ (2.18) 

which is an intrinsic statement written on the hypersurface S and the constant 
C depends on the spacetime curvature bound, only. 

Next, to estimate the gradient of u we consider the function 

V = v{Vu) := (1 + |Vup - fc)+, 

where k is chosen suitably large so that, thanks to the boundary gradient es- 
timate in Lemma 12.51 the function u(|Vm|) vanishes on the boundary of the 
hypersurface S. Multiplying (|2.18p by v'' for q > 1 integrating over the hyper- 
surface, and using Green's formula we obtain 

J (gt;9-i|Vwp-|-2u9|V2wp)dws 
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At this juncture, we observe that the higher-order term Au is controlled by 
the prescribed mean curvature equation (|2.16|) . We obtain 

J (^qv'^-^\Vv\^ + 2v'^\V^u\^^dvs<{q + l)C J dv^ 

and thus 

y(y(i+i)/^) ^ dv^ <{q + lfC f + w"-!) dv^. (2.19) 



To make use of (12.19^ it is convenient to return to our notation in coordi- 
nates, by observing that 

Vdetig) = Vl - Vdetig), 



so that 



|V(z;(?+l)/2)|2 > (lll)2„9-l|Vt;|2, 



where we used (if u > 0) 

C > ^l-\Vu\^ > C v-^'\ 

Therefore, provided we now assume that q > 3/2, (|27T9)) takes the following 
coordinate-dependent form: 

V(t;(?+^)/2)|2 dy<[q + ifC [ (««+3-l/2 + ^;9-l-l/2) dy . (2.20) 



By Sobolev's inequality in the local coordinates under consideration, we have 

/ r \ (n-l)/n /■ 

I u;2«/("-i)rfy] <c {\Vw\^+w^)dy, 

which we apply to the function w := ti('J+i/2)/2^ Recalling that r — 1 (after 
normalization) and observing that the domain of integration in y is bounded, 
we deduce from (|2.20p that 

/ r \ (n-l)/n 

[/ y(l+h)n/(n-l) dy) <C{q+lf («'?+3-V2+^,9-l-l/2)rfy 

or, equivalently, for all p > 2 

U dyj <C'/Pp^/P(j {v''+^+vP-^)dy] . (2.21) 
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Without loss of generality, we may assume that |lw||L°°(as) > Ij for otherwise 
the result is immediate. Then, (|2.21[) leads to the main estimate 

, / r \ ("-i)/(p") . 

It remains to iterate the above estimate, which yields 

\HL-in.)<C'\\v\\t^^^^^ (^^ vP«dy 



l/PO 



PO ^0 PO 

In conclusion, provided that po > 2n the sup norm of v is uniformly bounded 
by its L^" norm. 

Step 2. It remains to derive an estimate for some LP° norm. Following |14| . we 
return to the inequality (j2.16p satisfied by the function u and, for every A, we 
write 

A(e^") = Ae^"Au + A^e^^jVup 

> -CAe^"(z.(Vu))2 + A^e^^lVup. 

Combining this estimate with a direct calculation from (|2.18p (similar to the 
one in Step 1 above), we obtain 

A (t;«e^") >A2 e^^ + Aw^-^e^^ ((Vu, Vv) ~ Cv {v + 1)) 
+ gw«"ie^" (21 V^up - Cv^ + 2{Wu, V(Au))^ 



Then, by integrating over the hypersurface S, integrating by parts, using (I2.16P 
to control the term Au, and finally choosing A sufficiently large, we arrive at 



This completes the proof of Lemma 12.61 □ 



Proof of Lemma \2.7\ Step 1. We are going to control the sup norm of h, and 
to this end we will use Nash-Moser's iteration technique. Note that the elliptic 
equation satisfied by the second fundamental form a priori has solely L°° coeffi- 
cients. By scaling we can assume r — \. We consider an arbitrary point p' S S 
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and we set 5 :— d{p',dT,). Simons' identity [25] for the hypersurface reads 
Ahij = Ahij — {trh)ij 

= \h\^ % ~ {trh)hkhi,9^' ~ Tl^„ghklg^''g''' + R,p/,/i.fc5^«5'' (2-22) 
+ Vp(R,jN^)gP''- Vj(R,n), 

in which the Hessian {trh)ij vanishes since S has constant mean curvature. 
Recah here that N is the future-oriented normal to S. Thanks to (|2.22p we 
obtain 

A\h\' > 2 |V/.|2 + 2 \h\^ - C (|Rm|N + 1) \h\^ 
+ 2(Vp(R,JN^)5^''" V,(R,n),%). 

Let 93 be a smooth, non-negative, non-increasing cut-ofF function which 
equals 1 in the interval [0, 1/2] and in [1, 00). Then, the function -ijj := ip o k 
with K :— {d{p' , ■)/S) is a cut-off function on the CMC hypersurface S which 
vanishes near the boundary dT,. 

Fix some q G [l,oo). Multiplying p.23p by ip\h\'^ , integrating over E, and 
then integrating by parts, we arrive at 



> 



i^\h\'^A\h\Uvs 

(v|/i|'^(2|V/ip + 2 - C \h\^ - C |Rm|N(g + 1) |V/i|) 
-C\Rm\^\V^\\h\i+^)dv^. 



Using 



^Plhl"^ A\h\Uv^ < / 2\V7P\\h\'^+'\Vh\dvj:~ / 2qi;\h\'^\V\h\\'dv^ 
and Cauchy-Schwartz's inequality, we obtain 

<c£(^(^^+^^0«;|/,|«+2+(g+l)2^|/l|9 + i \^'ok\ l/j^+l^ dv^, 

(2.24) 

in which we can always choose ip so that < C |(y9|. 
Then, by Holder's inequality we have 

/ r \ 1/(9+4) 

ij^\h\i+Uv^] <CgS^-\ (2.25) 



In view of Lemma 12.61 the hypersurface is uniformly spacelike and so we have 
the Sobolev inequality 

" <C' [ |V(^|/ir+2)|d«s. 
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Combining this with (|2.24p and suitably choosing the function Lp, we find that 
for all i = 1, 2, . . . 



(2.26) 




Using the Nash-Moser's iteration technique we deduce that 



sup \h\<Cq6~^[ \h\'^^^dv^] . (2.27) 

B{p',S/2) \Jb{p\35/4) J 

Finally, choosing q^Q m (P?^ and g 2 in (P??7)l , we find 

C" 

sup \h\ < — . 

B{p',S/2} 

Step 2. Next, by relying on the sup norm estimate that we just established, 
we can estimate the covariant derivative of h. We need an estimate for the 
equation (|2.22p . From Gauss equation we see that the curvature of the hyper- 
surface is bounded by C'S~^ on the ball B{p',6/2). By introducing harmonic 
coordinates on the (Rieniannian) slice S, as in [18j. we see that the metric co- 
efficients belong to the Holder space C^'". Since the right-hand side of (|2.22p 
belongs to the Sobolev space W~^''^ for any q S (l,oo), thanks to the Sobolev 
regularity property for elliptic operators (in fixed coordinates) we find 

V^' C 




B{p',S/i) 



for some constant Cq > 0, which completes the proof of Lemma [^771 □ 



Proof of Lemma \2.8l We need now to estimate the time-derivative of the level 
set function u. Set p" = (u(7(sr)), 7(sr)) G S and consider the geodesic dis- 
tance function p = p(p",-) associated with the induced metric on the CMC 
hypersurface E. By the Gauss equation, the Ricci curvature is bounded (espe- 
cially from below) by 

Hence, thanks to the Laplacian comparison theorem, the distance function p is 
a supersolution for the operator —A + C" / p, that is, in the weak sense 

c 

Ap < —. 
P 
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Let ip be the (non-increasing) cut-off function introduced in the proof of 
Lemma [2771 Then, by differentiating with respect to t the equation (|2.12[) satis- 
fied by the solution u and in view of the bounds on h in Lemma 12.71 we obtain 

(a -\h\^- Ric(N, N)) (kV.) + . V'(^)) > 1 - 4 ^ 0' (2.28) 

where we have set e :— jC . Finally, applying the maximum principle (see 
(|2.30p below for the ellipticity property) we conclude that 

on = Bsr{^{sr), s^l'^rj^^ C {r sr}. 

□ 

Proof of Proposition \2.3\ and Theorem [KB. Step 1. The first variation CM{X) 
of the mean curvature along an arbitrary vector field X reads 

CM{X) ^ A{X, N) - {\h\^ + Ric(N, N)) {X, N) - {X, Vi?>, (2.29) 

where we recall that N is the unit normal vector field to the hypersurface and 
h is its second fundamental form. In the case of graphs, the linearization of the 
mean curvature equation around a constant mean curvature hypersurface reads 

£M{^) = A{iy{Wu) ifi) - + Ric(N, N)) i^{Wu) (p. 

Under our assumptions, this operator is uniformly invertible, since 

|Ric(N,N)| < C(7i)r-^ l/iP > > s-^/^r'^ (2.30) 

n 

and we choose s sufficiently small so that the term dominates Ric(N, N). 

More precisely, we have the following important conclusion. 

By the implicit function theorem, the linearized mean curvature operator 
CA4 on the space of all spacelike C^'" functions u, with fixed boundary value cr 
on dQs and with a € (0, 1), is locally invertible around any smooth hypersurface 
of constant mean curvature. In consequence, starting from any fixed spacelike 
hypersurface u with constant mean curvature 

2k{sr,r) G [k{sr,r), k{2s'^r,r)] 

and using the implicit function theorem, we can find a smooth family of spacelike 
hypersurfaces u* with constant mean curvature t varying in an e-neighborhood 
of 2k{sr,r). From Lemma 12.61 and Schauder's estimate, we then deduce higher- 
order uniform estimates for u*. So, the function u* is smooth and we can take a 
convergent subsequence of values t converging to the end-points of the interval. 

Next, by continuation we may still use the implicit function theorem and 
extend the smooth family under consideration for all mean curvature parameter 
values in the interval [k{sr,r), k{2s'^r,r)]. In other words, we conclude that 



.C'r'<^<-- 



dt 



C 
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for any t G [k{sr,r), k{2s'^r,r)] we can solve the Dirichlet problem (|2.12p of 
prescribed mean curvature equal to t, and, moreover, the solution u depends 
smoothly upon the mean curvature t. 

Next, differentiating with respect to t the conditions satisfied by the solution 
u we obtain 



A - {\h\^ + Ric(N, N))) (iy{\7u) ^) =1 in f7 



dt. 

Ou 

along d^s- 



(2.31) 



dt 

As already observed, by the maximum principle we have 

-Cr^ < 7^ < 0. 
ot 

This property shows that the family of CMC hypersurfaces forms a foliation of 
the region under consideration. 

Step 2. In this last part of the construction we choose the geodesic slice over 
which the CMC foliation should be based. We observe that the foliation con- 
structed in Step 1 need not pass through the given observer p = j{cr). To cope 
with this difficulty, we now vary the parameter s in order to ensure that the 
foliation contains a neighborhood of p. We proceed as follows. 

For any s G [c, 2c] we use the notation u^*^ for the function describing the 
CMC hypersurface constructed over the reference domain ilg C {t — sr} for 
the chosen value of the mean curvature 

2k{sr,r) € [k{sr,r),k{2s'^r,r)]. 

Now, we emphasize that u*^*^ depends continuously upon the parameter s. In- 
deed, in view of (j2.29p we can apply the implicit function theorem and we see 
that the solution depends smoothly upon the parameters arising in the domain 
of definition and upon the boundary values. Therefore, recalling the result in 
Lemma 12. 6[ given any function w^"^^ we may extend it to a whole family w^*) 
smoothly for all s G [c, 2c] . 

Then, by Lemma [2.4l we have u^"^ > T{s)r for some r(s) satisfying 2fc(sr, r) = 
k{T{s)r^r), which implies that t(2c) > c, at least when c is suitably small. 
Since we have u^''\'j{cr)) < cr, by continuity there is some so G [c, 2c] such that 
u(''°)(7(cr)) — cr. Therefore, we have constructed a family of CMC hypersurfaces 
u* with constant mean curvature 

t e [nk{sor,r),nk{2sQr,r)] for some Sq € [c, 2c] 

over some geodesic slice V.sg C {r = sqt}. Most importantly, the point jlcr) 
lies in the CMC hypersurface with mean curvature 2fc(sor, r). 
In addition, by a direct computation we obtain 

^ ' ^ \dtJ ^1 - |Vu|2 
and, in view of Lemma l2.8( the proof is now completed. □ 
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2.5 Further geometric estimates 

For any p' G St with Sr — d{p' , d'St) and thanks to our estimate of the second 
fundamental form in B{p',S/2) and Gauss equation, we see that the curvature 
of the hypersurface is bounded by C S^^r^^. By choosing harmonic coordinates 
as in [TB] and using the estimates in (|2.3ip . the function A = — i^(Vu) ^ 
satisfies (for any q G [l,oo)) 

\ 1/9 

^ ^ ■V^M'^dv^] <-^. (2.32) 



yB(p',5/4) J - 

In addition, let us investigate the geometry of the boundary dH of the folia- 
tion leaves. More precisely, we can estimate its second fundamental form IIoe, 
as follows. 

Proposition 2.9 (Boundary of the CMC foliation). The CMC hypersurfaces 
constructed in the proof of Theorem \2.^ also satisfy the uniform estimate 

r 

Proof. Recall that, for any tangent vector fields X, Y along dY,, the scalar 
IIas(X, y) is defined as g(Vx^, ^as), where Ngy^ is the normal vector field 
of 9S in the hypersurface E. 

On the other hand, since dH — dfl is obtained by the intersection of two 
level surfaces H :— {t = const.} and 5* := {d = const.}, dTi may be regarded 
as a hypersurface of codimension 1 in either Ti or S. The second fundamental 
form Ilg'j; of 91] in Ti. reads 

ll^AX,Y)=g{VxY,Nn) - i V?,d(X, F), (2.33) 

|Vd| 

where d is regarded as a function on Ti. and Vt^ denotes the covariant derivative 
associated with the induced metric on H, while N-h is the normal vector of dT, 
in Ti.. Similarly, we have 

lll^{X,Y)=g{VxY,N^) = J^V|f(X,y). (2.34) 

By a direct computation we find that 

V?,d(X, Y) = V'd{X, Y) - g(vd, |-) A{X, Y). 

So, we have 

^ ■g(x,y) < v2,d(x,r) < -^g(x,y) (2.35) 
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and similar inequalities for V|:t(X, F). On the other hand, by the triangle 
comparison theorem for the Riemannian metric g, there exists a constant C" 
(depending on c) such that 

This implies the existence of two functions a, b that are bounded by some uni- 
form constant C and satisfy Nqy^ = a N-u + h Ng. This completes the proof of 
Proposition El □ 



3 Local coordinates ensuring the optimal regu- 
larity 

3.1 Main statements for this section 

From now on we assume that the manifold M satisfies the Einstein vacuum 
equations. We will now prove: 

Theorem 3.1 (Local coordinates ensuring the optimal regularity). Given e > 
and q £ [1, oo) there exists a constant c{n, e) satisfying limg^o c(n, e,q) = Q such 
that the following property holds. Let (M, g, p, Tp) he a pointed Lorentzian 
manifold satisfying the curvature and injectivity radius bounds (|1.2p at some 
scale r > 0, together with Einstein field equation Ric = 0. Then, there exists 
a local coordinate chart x = (x") satisfying a;"(p) = 0, defined for all \x\'^ := 
(a:°)^ + (x^)"^ + . . . + (a;")^ < rf with ri :— ci(n, e)r, and such that 

sup (\gal3 - Vap\ + r \dgal3\) < e, 

f ^ \dh.p\'dx<C{e,q), 

' J\x\<ri 

where rj^^ is the Minkowski metric in these local coordinates. 

The proof of this theorem will be given at the end of this section, after 
establishing several preliminary results of independent interest. 

The main observation made in the present section is that the time function 
t associated with the CMC foliation constructed in the previous section admits 
well-controlled covariant derivatives up to third-order; cf. Proposition l3.4l below. 
Consequently, by following our arguments given earlier in [T^i. Proposition 9.1] 
we are led to the desired optimal regularity result in Theorem 13. II Recall that 
in the earlier work [12] we relied on a coordinate system in which the metric 
coefficients g had well-controlled first-order derivatives only; indeed, the time 
function in [12] was simply taken to be the geodesic distance function, which is 
only twice differentiable under the curvature and injectivity radius bounds. In 
contrast, in the present paper we have constructed a more regular time function 
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t (the mean curvature of the spacehke shoes) which turns out to have third-order 
regularity. 

Consider the constant mean curvature fohation St given by Theorem 12.21 
and observe that the time function t (together with the Lorentzian metric g) 
provides us with a natural flow <!>"' associated with the vector field 

-g(Vt,Vt)' 

such that the parameter m may differ from t by a constant. Denote by r the 
normal time function introduced in Section 2 and recall that, in the interior of 
the shce, 

a7<°' ^ -g(Vt,V.) '^-)<°' 

while the vector field Vi vanishes identically on the boundary. By starting from 
any arbitrary point pt £ 'St, the integral curve $™(pt) intersects each CMC slice 
exactly once and the flow $™ preserves the CMC foliation. 

Let y — (y*) be spatial coordinate chosen arbitrary on a given slice ^toi s-nd 
let us use the flow $™, in order to transport these coordinates to any other 
slice Et. Together with the mean curvature function t, these spatial coordinates 
provide us with spacetime coordinates y = {t,y^). The metric g then takes the 
form 

g = -X{t,yfdt^+g,,{t,y)dfdy^ (3.1) 
and satisfies the ADM equations 

2 A kij , 

(3.2) 

where A > is the lapse function and kij is the second fundamental form of 
expressed in the coordinates under consideration in this section 

From now on, without loss of generality we set r = 1. The central technical 
estimate of the present section concerns the lapse function and is stated in the 
following lemma. 

Lemma 3.2 (Second-order estimates for the lapse function). Under the as- 
sumption of Theorem \3.1\ and with the above notation, the function A satisfies 

For any S > we set Y,^ := {a: G S / d{x,dTj) > 6} and, away from the 
boundary of the slices, we can improve Lemma |3.2[ as follows. 



dt 



dt 
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Lemma 3.3 (Higher-order interior estimates for the lapse function). Under the 
assumption of Theorem \3.1\ and with the above notation, for any S > and 
q G [1, oo) one has 



sup 

£■5 



S-5 



(|fc| + |VA| + |V2A| 
(|Vfc|« + |v3A|« + 



dX 



dt 

V — 

dt 



dX 

dt 
d'^X 



d^X 



a<2 



<C{n,6), 



dvY, < C{n, q, S). 



Finally, based on Lemma 13.31 we prove that the time function t admits 
well-controlled third-order derivatives. Here, we use the covariant derivative 
V associated with the reference Riemannian metric in the coordinates under 
consideration, that is, g := X{t, yY dt^ -\- gij{t, y) dy^dy^ . 

Proposition 3.4 (Third-order estimates for the time- function). Under the as- 
sumption of Theorem \3.1\ and with the above notation, for all q G [l,oo) one 
has 

sup (\Vh\ + \Vh\ + I V'A|) < C{n, (5), 

J^^ {\VH\'^ + \Vh\^) dv^<Cin,q,6), 
where all the norms are computed with the reference metric g. 



3.2 Derivation of the key estimates on the lapse function 

This section is devoted to the proof of Lemma 13.21 

Step 1. Zero- order estimates in time. By integrating Weitzenbock identity 
(|2.17p and observing that |Ait| + |Vit| < C on a CMC slice E, we obtain 

/ iV^updvs < C+ / (V2?/)(nas,nas), (3.3) 

where uqy, '■= is the unit normal vector field of the boundary dT, on E. 
Since 

(V^M)(nas,nas) = Au - irg(V^-u) 

= An — |Vu| trg(Hgs) along E, 

we conclude with the boundary estimate in Proposition 12.91 that the second 
fundamental form k is uniformly bounded in the norm 

/ Ifcpdws < C. (3.5) 

Observe that this estimate covers the whole slice up to its boundary (in contrast 
with the interior sup- norm estimate given by Lemma l2.7p . 

We use the notation introduced in Section 2. Recall that the Riemannian 
distance function d = d{'j{s-\-s'^), •) takes the constant value cq = -f-s'^ on the 
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boundary 9S, and that cq — d is proportional to the intrinsic distance function 
to the boundary of the shce, i.e. 

<co-d< C"d(-,9S). (3.6) 

Moreover, by the Laplacian comparison lemma for distance functions and relying 
on our curvature assumption we have also 



C" ~ 1 

> Ad > on the hypersurface E. (3.7) 



Now, taking the trace of the second identity in (|3.2p and recalling that 
S = Et has constant mean curvature we obtain the elliptic equation satisfied by 
the lapse function 

AA = -1+ (|fcp + Ric(iV,iV)) A. (3.8) 
In view of (|3.7p and recalling that A > we deduce 

A(^\ + c^C'd^ > 

so that, thanks to the maximum principle, 

< A < c2C"(co -d). (3.9) 

In particular, in view of p.6p this implies the desired gradient estimate along 
tlie boundary at least 

sup|VA|<C". (3.10) 

as 

Next, by and the right-hand side of belongs to L^, which 

yields us a bound for the Laplacian of the lapse function 

\A\fdv^<C'. (3.11) 



On the other hand, by multiplying (j3.8p by the function A and integrating by 
parts, we find the gradient estimate 

/ iVApdws < C'. (3.12) 

Finally, by observing that 

/|AAp<C", AA|as=-l, i?ic(VA, VA) > -C" |VAp, 

integrating Bochner formula 

A|VAp = 2\y^X\^ + 2(VA, VAA) + 2i?ic(VA, VA), (3.13) 
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and then using (|3.10p together with a similar calculation as in (|3.4p . we arrive 
at an estimate of all second-order spatial derivatives of A: 

/ iV^Apd^s <C", (3.14) 

which is one of the estimates stated in Lemma 13.21 

Furthermore, we can also control certain nonlinear functions. Multiplying 
the identity 

AA^ = 2(|fcp + Ric(7V, TV)) A^ - 2 A + 2 |VAp, (3.15) 

by |VAp A"*^^"^^ on one hand and by |fcp on the other hand, for all e G (0, 1) 
we find 

/ A-i+'|VA|''dws < Ce, / l/cplVApdws < C. (3.16) 

Thus, by multiplying (|3.13p by |VAp and then using (|3.16[) . we obtain 

^ (iV^AplVAp + 1% V,Ap iVAp) dv^ < C, (3.17) 

where we used the Gauss equation for the expression of Ric. In particular, this 
provides us with a control of 



/ IVlVAI-^lrfws < C. 



Moreover, since the boundary values of |VA| are uniformly bounded, by Sobolev 
inequality we also have 

/ |VA|^ dws < C. 

Step 2. First-order estimates in time. This is the first instance where we use 
our assumption that the manifold is Ricci-flat. By differentiating (13. 8p in time, 
we obtain that 

KS) ^'^•') + + '"•=<'^''^') f + " 

- 2 ^ kkiKsg'^g'^g^^'X + Ric(V jlN, N) A + (V jLRic) (N, N) A 

at 3* ^ at ' 

+ t |VA|2 - 2//^fc,y VAVfcA + 2 A/' ViARfco, 

(3.18) 

where we used 

^flj = (v^Afcjj + VjAfc;, - V,A% + A(V,fc,j + Vjfc,, - V,%)) , 

= 2Ag"'5*^V,% + 2g'^ g'^'kuV ,\ - tg^'ViX 
= -2A5'='R,o + 2g''g'''kuV,X ~ tg^'ViX, 

(3.19) 
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as well as Codazzi equation Wikij — Vikij = Rizoj- Plugging in the vacuum 
Einstein equation Ric = 0, we arrive at the equation satisfied by the derivative 
of tie lapse function 



dX 



\dtJ 



dt 



= -Q + ViV', 

(3.20) 

in which, thanks to our estimates in Step 1, 

(iQP + IVI") dv^ <C, V'^ -2//^"%AVfeA. 

By multiplying (|3.20p by (ff )^^'^ on both sides, using Sobolev inequahty for 

the function (^)^^^^^j recalling ^ \dT.= 0, and finally integrating by parts, we 
conclude that 



C" 



< 



dt 



dt 



^2(2+^) 
2(l+e)\ 1/2 



,aA 

dt 



dX 



dt 



(n-2)/n 
I < 



1/2 







2 


d\ 




dX 


2 + € 














dt 




dt 




dt 





dvy 



1/2 



dX 



dt 



1/4 



\V\Uvy 



1/4 



(3.21) 

So, we should take e < ^-r if n > 5, but can take arbitrary e > if n < 4. We 
conclude that 

dX ^ 

dus < C if n > 5, 



dt 

dX 



dt 
dX 



sup 



at 



< Cg 

< C 



if n = 4, 
if n < 3. 



(3.22) 



In the case n < 3 above, we used once more Nash-Moser's iteration technique. 
Next, multiplying Bochner formula 



A ^9X 

by (co - d)^ using 



^2 9X 

dt 



.dX „ . dX 



2(Vij7,VA-) + 2i?zc V-,V- 



9t ' dt 



.dX _dX 



dt' dt 



,dX 
dt 



(co - d) 



dX 
'Tt 



dv^ < C, 
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and then integrating by parts, we find 



(CO - 



dvY. < C. 



(3.23) 



Multiplying p.lSp by |V^p and integrating by parts, we have 



<C + J |VA|' 
and, after using p.9p and p.23p . 



dX 

'at 




dt 


|VA| dvs 






at 


2 

dv-z 


+ 4 / A2 


— 

at 


2 



9i 



dvs < C. 



(3.24) 
(3.25) 



Now, multiplying Bochner formula (|3.13p by I^^P and then integrating by 



parts, we obtain 
(2 IV^Af 
< C 



ax 



at 



2i?ic(VA,VA) 



ax 



at 



'/JIVAP 



at 



W'x 



ax 



and, thanks to and ([X^ . 



IV^Al' 



ax 



at 



at 



h,v,x\' 



dvT, 

dvs + C' [ (AXf 



ax 



at 



dv^ 
(3.26) 



ax 



at 



dvT. < C. 



Step 3. Second-order estimates in time. We now turn to the most involved 
estimate concerning the function which, we claim, satisfies an equation of 
the form 

^( S) - IfcpS = V,F^ + /i + /2 + /3, (3.27) 



\at^ 



where 



at^ 

V'l^dv < C 



and /i G has the form /i ~ k*f[ (that a linear combination of such products) 
with 

' \f[\'<c, I 1/2!^ <C, 



and /a is bounded pointwise by C |V^Ap. This is one of the key observations 
in the present paper. 
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To establish (|3.27p we differentiate (I3.20p in time. It is not hard to show 
that all terms arising in the right-hand side of the equation, except those of 
form with \V\'^ < C, belongs to uniformly. We emphasize that we 

may arrange the other terms by introducing new terms V's so that they all 
have the desired form in ()3.27p . 

Let us now deduce from (I3.27P that 







1 










+ 


9*2 


1 



dws < C". 



(3.28) 



Namely, by multiplying both sides of (|3.27p by ^jy, then integrating by parts. 



and using |as= 0, we obtain 



< 



(l^l 9t2 

' a2A 



d^x 



9*2 



k\\f[\)dv^ 



(3.29) 



9*2 



I/2I- 



C 



92A 



9*2 



W^Xfdvj:. 



Since Ric > —Cg, by multiplying Bochner formula (|3.13p ) by |4t^I we obtain 



92A 



9*2 



.92A 



9*2 

U — 

9*2 



2(AA)2 
-f C|VA|2 



92 A 



9*2 
92A 



9*2 



|V|VAp| 
]WXAX\ 

By Cauchy-Schwartz inequality and Sobolev inequality, we then have 
2 

1/2!^ dvs 



dvj: 

dvY,- 
(3.30) 





_92A 


2 


92A 






+ 


9*2 



<C I \f[\''dv^ + C 
's 

' |VAAA|2 + |VA|4) 



C f |V2A|2|VA|2d«s (3.31) 



Hence, by combining with ([XS]) . (|3T7ll . (|3^ . (|3^ . and ([336]) we obtain 

(AA)4 



/ (|V2A|2|VA|= 



|VAAA|2 + IVAN dvs < C. 



which gives the desired estimate (|3.28p . 
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In summary, by combining the estimates already established in Lemma 12.71 
and in this proof, we thus have 



sup A |fc| 



dt 



|VA|^ + 1 dws 



+ 



|^(|V2Ap(|VA|^ 



dt 



V 



'dt 



dt 

d\ 

IVAI^ + IVAPIV— p 





dX 


2 


dkij 


2 




+ 


) 


dt 




dt 





dvs < C. 



(3.32) 



In the rest of this proof, we use the notation A Ki B when A ~ B is controled 
by the left-hand side of ((53^ . 
We can compute 



d_ 

dt 



^dX ,dX 

^dt I''! dt 

dX^ 
'dt) 
dX 
dt 



dX. . , , , , i9A ^ dkii , dX 



V, ( 2 A fc,, V, — ) - < ( V A, V — ) - 4 A fc,, fc.fe fcfc, — - 2 — (3-33) 



thanks to ()3.32p . and 

— (2 (RiNjNi kij)}?) — 4 (RiNjNi k^i)^~Q^ 8 RiNqNfcpgfcjpA'^ 



2( 



9t 



dt 



dt 



We can also compute the following: 

^(-4A(%,V,;V,A) 

dX dk ■ dX 

= -4 V.V, A) — - 4A V.V, A) ~ 4A , V.V, — ) 

- 8A fcy ViAVfeA + 4 A HVAp + 4A2 % V,/cfej V^A 
« Vi(4A^ %A:fcj VfeA) - SXkijkkj V^AV^A - 4A^ hjkkj V^VfeA, 



(3.34) 



(3.35) 



^^(2k,,kMkrs9''9'''g"'X^) 



dX 

4 kijkjkkki X + 12A kijkj^^kksks 



6 A TT—kjgksi 
dt 



0, 



(3.36) 
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(3.37) 

and finally 

^(-2//^%V,AVfeA 

= -8 A h,k,k V.AVfcA - 2^V,AVj A - 4fcy V,|^Vj A ^^'^^^ 
« V,(AV,|VAp) - AA|VAp « 0. 

To deal with the curvature term in (|3.33p we need some recall property of 
the curvature. We have 

with 

po („\ _ 1 po /„N _ 1 9A _ 1 %j 

^oo[S)-^Q^, ^o^[S)-^gy,, 2A2 9i ' 

T^fe \ ki 9^ T^fe ('„^ -'^ fci^5/i rk f\ j^k 

roo(g) -^5 rio(g) = -g — , rij(g) = r^, 

VjRifco; = VjRifco/ — r^j(g) Rofeo; ^ r"j,(g) Rioo; — r"Q(g) Hikai- 
Recall also the second Bianchi identity, 

VoRifcji — — VjRifcoi — VjRifcjo- 

We obtain 

— RiNjN = - V^^^'^J' j ■9''' + ^9^^ * * ^ 

= A,g^^V(RN***) + A.g^^ * * Rm + * VA * fc * Rn*** 
+ A Rn*n* * k * g~^. 



So, we have 

d_ 

dt 

which completes the proof of Lemma 



RiiVjAT, %)A^^ ~ 0, 



3.3 Proofs of the main statements 

Proof of Lemma \3.3[ Since the second fundamental form is bounded in each 
slice S*, then according to Gauss equation the intrinsic curvature of T,^ is also 
uniformly bounded by CS~^. Hence, from the injectivity radius theorem of 
Cheeger, Gromov, and Taylor for Riemannian manifolds [11], it follows that 
the injectivity radius of E*^ is uniformly bounded from below by C S. Next, 
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using the theorem of Jost and Karcher [TS], we can find a fixed number of 
harmonic coordinate charts covering S"^ and in which the metric is equivalent to 
the EucUdean metric and has W^''^ regularity for each q G [1, oo). In addition, by 
Sobolev's embedding theorem, the metric coefficients also belong to the Holder 
space C^'" for aU a e (0,1). 

Next, using an estimate from the equation (|2.22p satisfied by the second 
fundamental form in these harmonic coordinates, we deduce that kij £ VF^'? 
for all q. We also observe that the Christoffel symbols are of class C", so that 
this also provides us that Vfc S L*. All implied constants are uniform and only 
depend on the dimension n and the distance S to the boundary of the slice. 

Then, using a standard W^''^ regularity estimate for equation (|3.8p satis- 
fied by the lapse function (see, for instance, [16]) and noting that g € W^'^, 
we deduce that d'^X G L'^. Here, d'^X stands for any natural derivative in the 
constructed harmonic coordinates. Since dT G L'^, we obtain the third-order 
covariant derivatives |V'^A| G i^. Finally, we emphasize that for the spatial 
regularity of ^ and we need the crucial estimates established earlier in 
Lemma 13.21 We use an W'^''^ regularity estimate to equations (|3.20p and use 
Lemma [3.21 this leads to the desired estimate for Finally, the above esti- 
mates imply — |fcp^jT G W^^'"^, and we use again Lemma [3.21 and an 
regularity estimate in order to control □ 

Proof of Proposition \3.4\ We now want to control the covariant derivatives of 
the function t. Since this question is independent of the choice of coordinates, 
then on the fixed slice S'^ we choose finitely many spatially harmonic coordinates 
patches as in the previous proof, and we use them our new coordinates y*. Then, 
on this fixed time slice, the spatial metric belongs to W'^ '^ for all q G [l,oo). 
In particular, the Christoffel symbols r(g) — (Ff^) are uniformly bounded and 
dyr{g) G i'. Combining (|3.ip and Lemma 15^ together, we see that r(g) is 
bounded and |9t^j,F(g)| G L'^ at slice T,^ in these particular harmonic coordinates 
and at this fixed time. 

Now, we calculate the covariant derivatives of t in the coordinates chosen 
above. Since V^^i = — r(g)"^, we have sup-^s |V^i| < C. For the third-order 
derivative, we write V'^i — dt^yT{g) + r(g) * F(g) G L''. A direct computation 
yields us 

FOo(g) - FSJo(g) = 0, FO,(g) - FO,(g) = 0, F°.(g) - F°.(g) = ~^h„ 
~ dX 

Fgo = 2 A/-' -J , F^o (g) - Ffo (g) = 0, F^, (g) - F*:- (g) = 0. 
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Since V^t - V^i = -|fc, we have sup^s |V^t| < C5, and 

n 1 n 9A 1 9A 

VS.A = ^-if^V-A|VAP. 



A V 

Therefore, we have 



V i = (V - V)V^i + V{VH - VH) + VH 



and the result foUows from (|3.39p and Lemma [331 □ 

Proof of Theorem \3.1[ By a direct computation (see for instance [12 ) one can 
check that the Riemannian curvature of the metric g on IJ^ Sf is uniformly 
bounded and, actually, 

sup iRc^p^s - RcPjsI < Cs (IV^AI + \k\^ + |VA|2), 

hence sup^^^ 1^0/375! < C'^. According to the injectivity estimate for Riemannian 
manifolds established in [llj . the injectivity radius of the metric g at the point 
p is uniformly bounded from below, i.e. inj(M,g,p) > c^. Therefore, according 
to Jost and Karcher [TS] , we may choose harmonic coordinates a;" of g around 
P 

Noting that g = g + 2X^dt (g) dt, we obtain 

V^g =2 V^A^ (g) Vt (g) Vi + 2 A^ V^t (g) Vt 
+ 2 A^ Vi «) V^t + 4 A^ V^t (g) V^t. 

By combining with the result p.4p , we see that the coefficients of g in harmonic 
coordinates x" for the Riemannian metric g belongs to W'^'^ . □ 



4 CMC— harmonic coordinates of an observer 

4.1 Construction of local coordinates 
Preliminaries 

In Theorem l3.1[ we constructed coordinates in which the Lorcntzian metric coef- 
ficients have optimal regularity. However, one inconvenient of these coordinates 
is that they are not consistent with the CMC foliation constructed in Section 2. 
In the present section, we show that both strategies can be combined and we 
construct a new coordinate system which is based on the CMC foliation and has 
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optimal regularity as stated in Theorem 13.11 The basic idea is to now choose 
spatial harmonic coordinates on each CMC hypersurface. This strategy goes 
back to Anderson [1] who, however, used the time function given by a distance 
function. In contrast, in our present construction, the time function coincides 
with the mean curvature function of CMC slices and has much better regularity. 

In view of Theorem 12.21 since the second fundamental form is bounded in 
each slice Et and the spacetime curvature is bounded. Gauss equation implies 
that the intrinsic curvature of the slice is also bounded. So, according to [TT] 
there exists a constant rj = ri{n) > so that the injectivity radius of the slice 
is bounded below, that is, Inj(I](,pt) > 2rir, where pt is the orbit of the base 
point p along the above flow. By a theorem established by Jost and Karcher 
for Riemannian manifolds [18], there exists a constant rj' = rj'{n) > such that 
a harmonic coordinate system \y\ < rj'r exists around p on the slice '^t{p}: with 
p = (0, . . . , 0) and, on that slice, 

By using the above mentioned flow, the coordinate functions can be ex- 
tended to other slices Ej and, together with the time function t, yield a space- 
time coordinate system. Then, the Lorentzian metric g takes the form g = 
— A(f, y)^ dt^+gij{t, y) dy^dyK From the estimate of | Vtp given by Theorem l2.2l 
and in view of the expression v{S/u) ^ — —A, we deduce that 

Moreover, in view of the results in [18] and thanks to (|2.3ip and (|2.22|) . we have 
the uniform control 



Almost linear coordinates 

We now construct the coordinates of interest in this section. We can assume 
r = 1. For each i = 1, . . . ,n and for each shce Ef let a;' be the solution of the 
Dirichlet problem 

\x'=0 in E, n {y : |y| < r;'}, 
X = y on |y| = ?7 . 

Let t be such that the slice p G Ej has mean curvature t. 

By applying the maximum principle for the operator A we can derive some 
basic properties of the above functions. First of all, at the time t one has 

n 

n < At\y\^ = 2 ^ < 4n on the slice Ej, (4.4) 
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where we have solely used that the metric coefficients are uniformly bounded. 
Since 

— (^Atlyp^ = |VA*fc* VlypI < C5 on the slice Sj, 
and VA and k are uniformly bounded, we deduce from (|4.4p that 

I < Atlyp < 8n on any slice Ef n {|y| < 77'} 
for all |< - t\ < c,h"n - Fi'O™ 

now on we drop the subscript t in the notation. 
Then, by the same arguments as the ones above we find 

\Ay'\<e onS^nllyl <77'} (4.5) 

for all |t — i| < Now, since A(x' — y*) = — Ay% by the maximum principle 
we obtain 

C'er^'irj' - \y\) > - f > -C'er^'W - \y\) 

on the slice n {|y| < 77'} for all |t — t| < CsC"n • particular, along the 
boundary {|?;| = ry} the above property implies 

sup \V{x^ - 2/^)1 < Cin)ev' for \t - t\ < 
\y\=n' CsC'n 

Next, we can also estimate sup|j,|<^/ |V(x' — ?/')| from the equation satisfied 
by the coordinates, as follows. By integration by parts we obtain 

|V(x*-/)|2<C(n)e(77')"+^ (4.6) 

\y\<ri' 

Second, let w = max (O, |V(x* — y')P — C{n)eri'^ and consider Bochner formula 
A|V(x'-y^)P 

= 2|V2(x^ - y^)p + 2(V(.T^ - y% -VAy') + 2i?7c(V(x^ - z/), V(x* - y')), 



multiply it by w°- for a > 0, and integrate by parts. Then, using ()4.5p together 
with Sobolev inequality and Nash-Moser technique, we arrive at the sup-norm 
gradient estimate 



C f 

sup w < — ^ / wdy < C{n)ri'e 

\y\<n' ^ J\y\<v' 



(4.7) 



for all |t — t| < CsC"n ■ latter inequality follows from (|4.6p . 

By choosing e suitably small (depending upon the dimension only), (|4.7p 
implies that the harmonic map ^ = {x^ , . . . , a;") is a local diffeomorphism from 
{\y\ ^ v'} ^ onto its image. By the maximum principle, 5* is a map from 
{\y\ < v'} ^ to {|a;| < 77'}, which leaves invariant the boundary. Hence, "if 
is a diffeomorphism from {\y\ < 7;'} n Et to the Euclidean ball {|a;| < 77'}, and 
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) with {|x| < T]'} is a harmonic coordinate system. Moreover, by 
choosing e sufficiently small in (|4.7p we find 



1 1/9 9\ /9 9\ c) ( d d \ 



as inequalities between symmetric tensors. 
The ADM formulation 

Including also ^ t in the coordinates, we have therefore constructed local 
spacetime coordinates x^ , . . . , x") covering a neighborhood of the point p. 
Recall that N denotes the unit normal vector along slices St . 

Note that the function t appears as the time coordinate in two different 
coordinate systems, that is, ^' = {y°, y^, . . . , y") and ^' = (a;°, x^ , . . . , x") with 
y° = x° = t. It is easy to see that 

/ d \ _ d dx^ d _ d 



and so = A7V + ^, where we refer to ^ = X]"=i = as the 

shift vector. Define Aij — (f^j), and Ijki = gij{y,t)A~^ikA^^ji. It is not hard 
to see that dy^ — dx^ and dy^ — A^^ik (^dx'' — |p-rfa;"^ , hence in the coordinates 
(a;°, x^, . . . , x") the metric g has the form 



g = - A(t, y)2 dt^ + {t, y) dy'dy' 

= -X'^{dx")^+g,j{x",x){dx'+edx°){dx^ +^^dx"). 

For simplicity in the notation, we drop the tilde from gij and simply write the 
metric decomposition as 

g = -A2(da;°)2 + g^J{x°, x){dx' + Cdx°) {dx^ + f (4.9) 

Recall that the second fundamental form is defined by kij = (V a -^.N). 
where V is the covariant derivative associated with the metric g, and recall 
Gauss-Codazzi equations 

^ijki = Rfjki + kikkji — kilkkj, / , , „N 

(4.10) 

Vikij — Vikij — HuNj. 

The geometry of the slice is determined by the induced metric gij and the second 
fundamental form kij, both, satisfying the following evolution equations: 
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Note also that since x^, . . . are harmonic coordinates on Ef , we have 

92 



where Qij{dg, dg) is some quadratic expression in dg with coefficients depending 
on the inverse metric g~^. 

Estimating the shift vector 

Next, we derive the equation for the shift vector ^. By differentiating the har- 
monic equation As*^ = with respect to x", and using (|4.11[) . we get 

= g'^'g'^ (v,(-2Afc;j + + ViQ + V,{-2\ku + + V^^,) 



= 2 {^i^ + g^'R,j^^ + g''^Vi{XtrK) - 2g^^ g'^ ku\I jX - 2X{trk)i + 2Xg''^RiN 

where A^'^ is the fc-th component of A^. By combining this result with the 
constant mean curvature equation, this gives us the elliptic equation satisfied 
by the shift vector 

AC" = -g'^R^je - {trk)g''' WiX + 2g''' g'^ kuW , X - 2Ag*^'Rzw. (4.13) 

It is easy to see A|^| > —C for some constant C depending only on the dimen- 
sion. By choosing a sufhciently large constant we obtain A(|^| -I- C" |a;p) > 0, 
hence by the maximum principle we arrive at the following sup norm estimate 
for the shift vector 

\^\<Cin)W~\x\). (4.14) 
4.2 Proof of the main theorem 

We are now in a position to give the proof of Theorem 11.11 By scaling, we may 
assume r — 1. 

Step 1. Spatial derivative estimate. We are going to use (|2.22p (|4.12p . (|3.8p . 
(|4.13p . together with elliptic regularity estimates, and establish a bound for the 
spatial derivatives of the metric g. 

By choosing other harmonic coordinates on each slice, letting 77' suitably 
small, and recalling the U' regularity estimates for uniformly elliptic operators, 
we find for all q G [1, 00) 



/ \V\x^-y%''<C,. 

J\x\<ri' 



(4.15) 



This implies that 

< Cq (4.16) 



/ 


dgtj 


'\x\<7j' 


dx^ 
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and, therefore, for all a £ (0,1), |lg|lc°{|2;|<?7'} < Cq- In view of (|4.12p . we 
obtain 



x\<r]' 



9 



< c„ 



(4.17) 



and using estimate, since the coefficient of the Laplacian operator are Holder 
continuous and the function under consideration vanishes on the boundary. 



d 



< 



C{n) 



X I 



92 



dx'^dx 



((^"-N').9., 



(4.18) 



Note that in the expression = 5''' g^kg^i + T * + cT * ^, we have T e L« 
(thanks to (|4T6)) ) and 

\drH\<C{\d'g\ + \dg\'){rj''^\x\') 

thanks to (|4.14p . The latter term belongs to in view of (|4.18p and therefore 
LP regularity estimates applied to (|4.13p yield 



SUp|9a;^| 



\die\'<c. 



x\<rj' 



or, in covariant form, we have estimated the first- and second-order derivatives 
of the shift vector 

(4.19) 



sup|vei+/ iv2er<c, 

\x\<r]' 



In addition, since dxk = Vfc + T * k and |Vfc| € L'^ by Lemma [3.31 we also 



find 



Similarly, since d^X ^ V^A + F * VA and in view of Lemma [331 we also obtain 

x\<ri' 

In summary, we have now control the spatial derivatives (up to second order) 
of the metric, the lapse function, and the shift vector: 



{v''~\xr)g^„X,e&W'^''{{\x\<lf}). 



(4.20) 



Step 2. Estimates of first- order time derivatives. The strategy now is to differ- 
entiate the equations (I2.22p , (|4.12p , (|3.8p , and (|4.13p with respect to t and then 
use the elliptic regularity property. 
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First of all, thanks to Step 1 we have C^gij — € W^''', A kij S 

W^'i for aU q G [l,oo). By we have ^ G W^^'' and , in particular, 



dxdx^ 



^ G L|, i.e. in other words for all g G [1, oo) 



sup 

|a:|<i)' 







d^gi] 




•J \x\<ri' 


dxdx^ 



< Ca 



(4.21) 



In view of Step 1 and (|4.1ip again, we have G £| for all q G [l,oo). 
Then, from Lemma 13.31 we deduce 



sup 

\x\<ri' 



dX 



dt 



Since - 



dt 



9t2 



|Vfc|^ 



+ |V'A| 
dX 



V 



X I 

g^ = § + (e,VA), wehave 



dt 



V 



(4.22) 



and 



„ 9A ^dX , 

= V— + Ve * VA + e * V^A 



f9A 8X 

= V^— + V^^ * VA + C * V^A + * V^A. 
ax'^ ot 



Then by gH]), and Step 1, we find 

, i9A , dX . 



\x\<r)' 



(4.23) 



Next, note that 



dg 



By differentiating (|4.13p with respect to the time variable we obtain 



A 



A1 + A2 + A3, 



(4.24) 



with 

A2 := 



25 5^a^V,A, 



a D 

- 2A,-(V^Ric)(^,A.) - 2A,-Ric(V^^,A^) 



A 



2Ag'^'Ric(^, V_^7V) + ^ (g'^-g-i?,^^^- + trfc /V' V, 
25'=''/'.g^^fci,Vj A - 2g'='g^'-g^"^fci,VjA + 2A5'=''g'^RzAr + g''' g''^ 3^] 
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and 

Since the spacetime under consideration satisfies the vacuum Einstein equa- 
tion, we obtain 



dRic T ^ dg . 9 



+ ^ * (^g~'^ *Ric*^ + k*g~^*S/X + 5"^ * V^^) . 



(4.25) 



It is a classical observation that 



-2S = A.('fe)+V,y,+V,F„ (4.26) 



dxo \ ax" 

where 

is the Lichnerowicz Laplacian, and 



2 aa-o^ y ''aa;0" 

Since e Lg, = + T * V^, we see that ^ G W-i'" for all 

g e [1,00). Note that ^ e W^.-z, |^ g i|,VA e V^V.^'^ e L«,V^ e 

i^, V(^) G L|, V2(^) = 5V(^) + r * V(^) e 
Now at the boundary \x\ = 77', we have 

where wc used ^ 0. By applying the estimate, we conclude that 

§1^ G Vl/jt''^ for all q G [l,oo). In particular, 

= V^g; + V(0) . p-> . V« + p-> . V^(g,) . « . H'-^.. (4.27) 

In summary, we have proved that the first order (in time) derivatives of the 
metric, lapse function, and shift vector have well-controlled spatial derivatives 
up to first (or even second) order: 

IfeH--', §^,W}'. '£,Wi' ,4,28) 

for all q G [1, 00). 
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Step 3. Second- order time derivative of the metric and lapse function. 
First of all, by differentiating (|4.1ip we find 

d'^h, .jdX ^ dk , „ f d£, 

dx° dx" ^dx° dx^ \dx" 

dk _^ _2 dg 2 dRic 



(4.29) 

Recalling that g,,,e, A € W^^,,^ ^ g ^ g ^2,,^ ^ ^ ^1..,^ ^ ^ ^1^,^ 

^ e W-^''i, and combining together with gTH), g^S), and (H^Ol), 
we get the following bounds for the metric and the second fundamental form 



dx"^' dx° 

for all q E [I, 00) 
To handle tl 
XN + ^, so that 



^eL«, ^ew-'-^ (4.30) 



To handle the lapse function we note that — §F + (Cj^'^) and 



dx°^ dx"^ ' dxO 
^(AiVA) = ^+e(^)-e(e,VA) 

a^A „ dX 



dt^ ' (e^V— )-(V,e,VA)-V^A(e,0, 



and 

a^ 



9^^ * * C * VA + .g"^ * ^ * * VA + 5"^ * ^ * ^ * V^A 

OT*^ ■ OX^ ox" 

g^^ * ^ * V^A + g^^ * * VA + g^^ * V^A * 

ox" ox" 

-2 ,,, V72(- ,,, C ,,. YT \ I „-2 ,,, V7 e ,,, V7 e ,,, V7 \ I „-2 ,,, Y7(- ,,, f Y72 i 



g^^ * V^C * C * VA + g^^ * * * VA + * V$ * ^ * V^A 

(4.32) 



1 ^> ^ (^-^ 1 , 5A 
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Hence, combining together (|4.19p . (|4.22p . and (|4.28p . we arrive at the fol- 
lowing control of the lapse function 











/ 


V 9 


+ 




J \ x\<.7]' 







for all g e [1, oo). 

Step 4- Second-order time derivative of the lapse function. 

It remains to derive the second-order time estimate for the shift function. 
By differentiating (I4.24p in time, we have 



(4.33) 



with 



d^S,^ ( dS, dg _2 , _3 ,dg ^ \ 



-2 , 

+ q * k* V 



5.g 



+ 9 * 



-3 ^9 dk 



dg 



:= 



d^Ric 



1 ^ dRic f d£ 1 n do _ 



dxo 



dg ,2 , % , 9g , „2_^ 



+ 9-''V'(Sr).S7 + V(^).!7-='.V 



*9^ *^*g~^ 
dx° 



Note that we already have dx,tg, jg G for all q E [1, 00), except that we 
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a2c 

do not control yet. Therefore, we can write 



A = — g Hii n+g * ^ * VA 



^9:e0 ^ dx^^ ^dx^' dx° 

g-^*g* V^f-^) * e + V(-^) * * Ve mod. 



In view of g2Ql), (|4:28l) . ([4271 . and (jOO]) we have 

dx^ ' dx°^ dx°^ dx^ ^dx°^ dx^' 

dx^^ 



(4.34) 



and 



g'^ * fv^ — ^ +Rni* 



aa;0^ V dxO' dx°' 



for aU g G [1, cx)). Consequently, we have 

dxO^ ^dx°^ « n^m s ; ^.J J ^^ gg^ 



where for fixed k, etc. are tensors, and V and d are covariant derivatives 
and partial derivatives in the coordinates x% respectively, with moreover 



x\ <?7' 



For the second equality, we used S^d^g G L'^. Now, we will use the L'p regularity 
estimates in the following manner. 

Since the coefficients of the elliptic operator g'^''dadb belong to C*" on the 
closed ball {\x\ < 77'}, we can solve the equation g"-^ dadbu'^^^ = F^^^ on \x\ < rj' 
with the trivial boundary condition u'^^^ \\x\=rt'— 0. We then apply the 
regularity estimate and obtain 



\x\<rj 
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Next, we observe that 

+ d{dg * d'^u^,^^ * ^ + a^yfcjJ * a^) + d^u'^^ * * g-^ 
and so, for some new terms F*"'' , we obtain 

+ = .g°'a„9b(7.^r ) + d.F'^' + (4-36) 

for some 



IV < C,. (4.37) 



x\<.ri' 

Since 

= g'^'dadtiutD + d{r*u*o + d^9*nH, 

and £,d^g S i'' for all q, by modifying i^''* and F*^, we can show 

A (0J - z.f„r ) + (£1 - <r ) - + F^ (4.38) 

where the notation A stands here for the covariant Laplacian of a vector field. 

It is not hard to see \\x\=-q'— 0. Let w"' = — Integrating 

^|^fe|25+2 lY^Q induced (intrinsic) volume form, using v*' \\x\^ri'= a-nd 

(|4.38p . and finally applying Holder inequality, we find 

/ |f'=p9|Vi;'=p < C, / |«'=p9+2 + C, (4.39) 

for all q G [1,cxd). This implies, in particular, /|^|<;^/ < C by Sobolev 

inequalities. Combining this result with (|4.37p . we arrive at the estimate for the 
shift vector 

< 

In summary, we have obtain the following uniform control of the second- 
order time derivatives of the metric, lapse function, and shift vector: 

^^Ll, ^eLl ^eLl (4.40) 
dx^^ dx°^ 

for all q G [l,oo). By combining ((i?^ with and ((i^ . the proof of 

Theorem II. II is now completed. 





Q2^k 


J \ x\ <7]' 


dxo^ 



45 



Acknowledgements 



The authors thank Lars Andcrsson for providing them with bibhographical 
informations. The first author (BLC) was partially supported by Sun Yat- 
Sen University (Guangzhou) through a grant "New Century Exeellent Talents" 
(NCET-05-0717). The second author (PLF) was partially supported by the 
Centre National de la Reeherche Scientifique (CNRS) and the Agence Nationale 
de la Recherche (ANR) through the grant 06-2-134423 entitled "Mathematical 
Methods in General Relativity" (MATH-GR). 



References 

[1] M.T. Anderson, On stationary vacuum solutions to the Einstein equations, 
Ann. Henri Poincare 1 (2000), 977-994. 

[2] M.T. Anderson, On the structure of solutions to the static vacuum Einstein 
equations, Ann. Henri Poincare 1 (2000), 995-1042. 

[3] M.T. Anderson, On long-time evolution in general relativity and geometrization 
of 3-manifolds, Commun. Math. Phys. 222 (2001), 533-567. 

[4] M.T. Anderson, Regularity for Lorentz metrics under curvature bounds. Jour. 
Math. Phys. 44 (2003), 2994-3012. 

[5] L. Andersson, Constant mean curvature foliations of flat space-times. Comm. 
Anal. Geom. 10 (2002), 1125-1150. 

[6] L. Andersson, Bel-Robinson energy and constant mean-curvature foliations, 
Ann. H. Poincare 5 (2004), 235-244. 

[7] L. Andersson and V. Moncrief, Elliptic-hyperbolic systems and the Einstein 
equations, Ann. Inst. Henri Poincare 4 (2003), 1-34. 

[8] L. Andersson and V. Moncrief, Future complete vacuum spacetimes, in 
"The Einstein equations and the large scale behavior of gravitational fields", 
Birkhauser, Basel, 2004, pp. 299-330. 

[9] R. Bartnik and L. Simon, Spacclikc hyporsurfacos with proscribed boundary 
values and mean curvature, Commun. Math. Phys. 87 (1982), 131-152. 

[10] A. Besse, Einstein manifolds, Ergebenisse Math. Series 3, Springer Verlag, 1987. 

[11] J. Cheeger, M. Gromov, and M. Taylor, Finite propagation speed, kernel 
estimates for functions of the Laplace operator, and the geometry of complete 
Riemannian manifolds, J. DifF. Geom. 17 (1982) 15-53. 

[12] B.-L. Chen and P.G. LeFlocii, Injcctivity radius estimates for Lorentzian 
manifolds, Commun. Math. Phys. 278 (2008), 679 713. 

[13] D.M. DeTurck and J.L. Kazdan, Some regularity theorems in Riemannian 
geometry Ann. Sci. Ecole Norm. Sup. 14 (1981), 249-260. 



46 



[14] C. Gerhardt, H-surfaces in Lorentzian manifolds, Commun. Math. Phys. 89 
(1983), 523-533. 

[15] C. Gerhardt, Curvature problems, in "Series in Geometry and Topology", vol. 
39, International Press, Somerville, MA 2006. 

[16] D. GiLBARG AND N.S. Trudinger, Elliptic partial differential equations of sec- 
ond order, Springer- Verlag, Berlin, 1983. 

[17] S. Hawking and G.F. Ellis, The large scale structure of spacetime, Cambridge 
Univ. Press, 1973. 

[18] J. JOST AND H. Karcher, Geometrische Methoden zur Gewinnung von a-priori- 
Schranken fur harmonische Abbildungen, Manuscripta Math. 40 (1982), 27-77. 

[19] S. Klainerman and I. Rodnianski, Ricci defects of microlocalized Einstein 
metrics, J. Hyperbolic Differ. Equa. 1 (2004), 85-113. 

[20] S. Klainerman and I. Rodnianski, Rough solutions of the Einstein- vacuum 
equations, Ann. of Math. 161 (2005), 1143-1193. 

[21] S. Klainerman and I. Rodnianski, On the radius of injectivity of null hyper- 
surfaces, J. Amer. Math. Soc. 21 (2008), 775-795. 

[22] S. Klainerman and I. Rodnianski, On the breakdown criterion in general 
relativity, preprint, 2008. 

[23] B. O'Neill, Semi-Riemannian geometry with applications to relativity, Acad. 
Press, New York, 1983. 

[24] R. Penrose, Techniques of differential topology in relativity, CBMS-NSF Region. 
Conf. Series Apph. Math., VoL 7, 1972. 

[25] P. Petersen, Convergence theorems in Riemannian geometry, in "Comparison 
Geometry" (Berkeley, CA, 1992-93), MSRI Publ. 30, Cambridge Univ. Press, 
1997, pp. 167-202. 

[26] M. Reiris, The constant mean curvature Einstein flow and the Bel- Robinson 
energy. Preprint larXiv:0705.3070] 

[27] M. Reiris, The ground state and the long-time evolution in the CMC Einstein 
flow. Preprint T M"Xiv:0809.3444] 

[28] J. Simons, Minimal varieties in Riemannian manifolds, Ann. of Math. 88 (1968), 
62-105. 



47 



